Spin bottleneck in resonant tunneling through double quantum dots 
with different Zeeman splittings 
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We investigated the electron transport property of the InGaAs/GaAs double quantum dots, the 
electron g-factors of which are different from each other. We found that in a magnetic field, the 
resonant tunneUng is suppressed even if one of the Zeeman sublevels is aligned. This is because the 
other misaligned Zeeman sublevels limit the total current. A finite broadening of the misaligned 
sublevel partially relieves this bottleneck effect, and the maximum current is reached when interdot 
detuning is half the Zeeman energy difference. 

PACS numbers: 73.63.Kv, 72.25.Mk, 73.23.Hk 



Electron g-factors in III-V semiconductor heterostruc- 
tures can be tuned by changing the alloy ratio and thick- 
ness of each quantum well. Novel spin-related physics 
in such a g-factor-tuned system has attracted consider- 
able interest in the past decade [1]. Electrical tunings of 
electron g-factors has been demonstrated in single and 
coupled double quantum wells. The g-factor tunings are 
accomplished by changing the position of the electron 
wave functions that is spatially delocalized over the re- 
gions with different g-factors [2|,[3,y]. The g-factor tuning 
used above is a powerful tool for the coherent manipu- 
lation of the electron spins, however these works were 
done for an ensemble of spins in quantum wells. To re- 
alize the spin based quantum information devices, which 
is a one of the most ambitious applications of semicon- 
ductor spintronics, the coherent manipulation in a single 
spin level is necessary [5]. The quantum dot is known 
as a system where a number of electron as well as their 
spin state can be well defined and easily controlled [6|. 
An application of such g-factor tuning to quantum dot 
systems offers a novel candidate for future electric de- 
vices for coherent spin manipulation. Tuning the Zee- 
man splitting gifisB of each spin in a quantum dot ar- 
ray has been proposed for individual addressing of spin 
qubits and fast gate operation between two qubits, where 
gi {i = 1,2,---) is different for each dot [HI, Q- Selec- 
tive addressing has been demonstrated so far in a double 
quantum dot made from spatially homogeneous g-factors 
with an additional micro-ferromagnet nearby. The mag- 
net creates a ^lOmT field difference of external magnetic 
field B in each dot [1]. Quantum dots with different 
can offer a much larger Zeeman energy difference in the 
same external field with smaller spatial separation. 

In this letter, we report a novel behavior of a single spin 
in a g-factor-tuned double dot, where Zeeman splittings 
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differ greatly from each other in a magnetic field. We 
investigate this system in a simple, well-defined regime 
where the total electron number N in the double dot is 
zero or one. Accompanying the theoretical calculation, 
we reveal that the resonant tunneling via two Zeeman- 
split levels is suppressed even if one pair of Zeeman-split 
levels is aligned. This novel spin bottleneck effect is par- 
tially relieved by a finite broadening of zero-dimensional 
states owing to interdot and/or dot-electrode tunnel cou- 
plings. As a result of competition between the bottleneck 
and the level broadening effect, the current is maximum 
in the configuration where interdot level detuning is set 
to half the Zeeman energy difference. 

Vertical double quantum dots with different g-factors 
are formed in a submicron-scale pillar of a triple bar- 
rier structure with a surrounding Ti/Au Schottky gate, 
as schematically shown in left inset of Fig. 1. The 
structure consists of seven layer structure including 
triple barriers and comprises, from top to bottom (or 
from left to right in the schematic potential energy 
landscape in the right inset of Fig. 1), a gradi- 
ently n-doped Alo.05Gao.95As source electrode, 7-nm- 
thick Alo.3oGao.7oAs barrier, 7.5-nm-thick Ino.04Gao.96 As 
well, 6.5-nm-thick Alo.30Gao.70As center barrier, 10-nm- 
thick GaAs well, 7-nm-thick Alo.30Gao.70As barrier, and 
gradiently n-doped Alo.05Gao.95As drain electrode. The 
fabrication prosedure is the same with the previous ver- 
tical double dots [9] . Measurements were performed in a 
dilution refrigerator at an effective electron temperature 
of ~0.1 K and in magnetic fields of up to 15 T applied 
perpendicular to the wells [10]. 

Figure 1 shows the differential conductance, 
dlsD/dVsD^ plotted as a function of source-drain 
voltage VsD and gate voltage Vq in zero magnetic field 
pj| . Current steps are recognized as dark blue lines. 
In the positive Vsd region, near the current threshold, 
several current peaks (not steps) appear, as marked by 
arrows. These peaks are due to the resonant tunneling 
through the ground state of the left dot and an excited 
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FIG. 1: dIsD/dVsD plot of the sample measured at T. The 
arrows mark current peaks due to resonant tunneling where 
the ground state of left dots is aligned to one of the excited 
states of the right dot, as shown in the right inset. Left inset: 
schematic of a gated vertical double quantum dot structure 
with different g- factors for each dot. The two quantum dots 
are made of InGaAs and of GaAs layers. We apply Vsd to 
the substrate-side (bottom) electrode and measure the drain 
current from the top electrode. Right inset: potential energy 
landscape for the resonant tunneling condition. 



state of the right dot, as shown in the right inset. These 
current peak lines run nearly parallel to the Vq axis 
because the side gate capacitively couples to the two 
dots almost equally, and the alignment of the two dot 
levels is nearly maintained against Vq- These behaviors 
have also been observed in vertical double dots with the 
same quantum well layers f^. In the first Coulomb 
staircase defined as the TV = threshold line and the 
neighboring parallel line, the current is carried by three 
charge configuration {Ni^ N2) = (0,0), (1,0), and (0,1), 
where Ni {i = 1,2), is the number of electron in each 
dot. In this configuration, the situation is simple enough 
to avoid electron-electron interactions in each dot. 

Application of the perpendicular magnetic field change 
the energies of excited orbital states of the dots, and shift 
Vsd at the resonant tunneling peak lines [l2|. Hereafter, 
we concentrate on the current peak line found around 
Vsd "^35 mV and carry out detailed measurements un- 
der various magnetic fields. Figures 2(a)-(d) show the 
diDS /dVsD plots for several magnetic fields. The current 
peak line, which is recognized as adjacent blue and red 
lines, have a clear kink structure. In vertical double dots 
with the same quantum well layers, current peak lines 
are always straight in any magnetic field and the kink 
structure has never been observed. The kink is charac- 
terized by two values. Si and J2, as marked in Fig. 2(d). 
Both Si and S2 linearly increase with magnetic field, as 
shown in Figs. 2 (e) and (f). With low magnetic fields 
the peak line becomes straight asymptotically. Similar 
kink structures are found in all other current peak lines 
in various positive Ksd's at high magnetic fields, as ling 
as that the width of the peak is narrow enough to resolve 
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FIG. 2: (a)-(d) Differential conductance, dlsu/dVsD, as a 
function of source-drain voltage and gate voltage under sev- 
eral different magnetic fields, showing increasing kink struc- 
ture characterized by ^1 and S2. Magnetic field dependences 
of (e) Si and (f)fc. 
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FIG. 3: Schematic diagram of resonant tunneling peak line 
for Zeeman mismatched system, and characteristic potential 
landscapes A-D. 



kinks. 

In a double quantum dots with the same g-factors, 
both Zeeman sublevels are aligned at the same time at 
a certain Vsd^ regardless of the magnetic field [12]. In 
double quantum dots with different g-factors, however, 
the alignment of each Zeeman sublevel is achieved at 
its own VsD^ as schematically shown in Fig. 3, in ac- 
cordance with the conditions of the measurements. We 
assign four particular conditions labeled A to D in the 
magnetic field. Under condition A, the aligned levels for 
up-spin for both dots is in the transport window and 
the resonant tunneling current is carried out by up-spin 
electrons. By increasing Vg^ i-e.^ by lowering the energy 
level of both dots, the down-spin Zeeman sublevel in the 
left dot comes within the transport window under condi- 
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FIG. 4: Calculated current by the Bloch equation method, 
(a) Differential conductance plotted as a function of inter- 
dot detuning (normalized by Zeeman energy difference A = 
\ql — qrIiubB) and left-dot potential (normalized by left-dot 
Zeeman splitting energy IqlIij^bB), where 2. 5 A = lOQ = 
IOOO7L = IOO7H = 7. Transport condition is categorized 
as two cases, I and II. We neglect a finite transition region 
between I and II due to small here.(b) Peak current in 
cases I and II as a function of normalized interdot coupling 
^1/7. Other parameters are the same as in (a). 



tion B. Under this condition, although even the up-spin 
states are aligned and the resonant tunneling channel ex- 
ists, the resonant tunneling is Coulomb blocked once the 
down-spin state in the left dot is occupied. We name 
this suppression process the spin bottleneck. Similarly, 
under condition C, the occupation of the up-spin in the 
left dot prohibits the subsequent tunneling even though 
the aligned down-spin channel exists. Thus under both 
conditions B and C, the bottleneck channel is occupied 
and there is no steady-state resonant current. It should 
be stressed that this bottleneck effect is similar to the 
Pauli spin blockade [13, 14, 15], since the stochastic single 
electron occupation of the bottleneck channel ultimately 
leads to a blockade in both cases. 

The bottleneck can be lifted by level broadening, which 
is induced by finite tunnel couplings among the dots 
and the electrodes. The broadening couples the mis- 
aligned Zeeman sublevels and provides an escape path 
for the electron in the bottleneck channel, and relieves 
the bottleneck effect. The electron transport is carried 
out within the competition between the bottleneck and 
the escape effects. While smaller detuning in the reso- 
nant tunneling channel can carry large current, smaller 
detuning of the bottleneck channel is required to avoid 
the bottleneck. As a result of the compromise under the 
intermediate detuning condition D, not under B or C, 
the maximum current is expected. 

This scenario is supported by the theoretical analysis. 
We evaluated the resonant tunneling current through two 
dots with different Zeeman splittings by the Bloch equa- 
tion method [H, [l^- The quantum dot states are de- 
scribed by five bases, |0,0),| T,0),| j, 0), |0, T), and |0, j). 
Each dot has own g-factor, and gu and each level has a 
Zeeman sublevel and their difference isA= IgL—gnlfJ^BB 
(hereafter, we assume Ql < 9r < 0). In the calculation, 
we neglect the cotunneling and tunneling processes with 
phonon absorption/emission and assume zero tempera- 



ture. We fix VsD to be large enough so that the right 
Fermi level is far below the right-dot level cr. Figure 
4(a) shows the differential conductance plotted as a func- 
tion of cl — and cl- In the figure we set the origin 
of cl — ^R at the position where the lower Zeeman (up- 
spin) sublevels are aligned and is measured from the 
Fermi level of the left electrode. There are two conditions 
labeled case I and case II. In case I, only the up-spin 
Zeeman sublevel of the left dot is within the transport 
window (corresponds to scheme A in Fig. 3). In case 
II, both the up-spin and down-spin levels of the left dot 
are within the transport window (schemes B-D in Fig. 
3). In both cases, we obtain a current with a single peak 
of Lorentzian shape. In case I, the current has a peak 
at zero interdot detuning. In case II, the maximum cur- 
rent is shifted in the negative detuning direction by the 
amount of A/2. These peak positions in both cases I 
and II are found to be independent of parameters such 
as the coupling of the left (right) dot to its neighbor- 
ing electrode, Jl{r)^ and interdot coupling Q. Note that 

— and cl roughly correspond to Vsd and Vg in 
the measurement and Si and 62 correspond to \gL\fJ^BB 
and A/2, respectively. Thus the clear kink structure in 
the current peak line in Fig. 4(a) is similar to those in 
Figs. 2(a)-(d). Peak currents are plotted in Fig. 4(b) as 
a function of Q. The peak current in case II vanishes for 
small Q. In contrast, the two peak currents take simi- 
lar values for large where the dot-electrode couplings 
limit the current. In Figs. 2(a)-(d), the current peak 
heights are almost the same (~10 pA), both below and 
above the kink structure. This corresponds to the large 

Indeed, our previous studies on vertical double dots 
with similar barrier thicknesses indicate Q ~0.1 meV and 
1l{r) ^ O.OlmeV [9]. In this sample, 7^ is much smaller 
than jR [11]. In Fig. 2 (a)-(d), transition between the re- 
gions corresponds to the case I and II is not as abrupt as 
in Fig. 4 (a) due to the finite temperature. It should be 
noted that, although this calculation reproduces the shift 
of the peak and the peak height well, the calculated peak 
width is always larger than the peak shift A/2, whereas 
in the measurement, A/2 seems to be twice as large as 
the peak width at 13 T. This discrepancy might be solved 
if we include the cotunneling effect. 

Here we evaluate g-factor of each dot. In order to con- 
vert from 5i and 62 to \gL\fiBB and {\gL - gR\)fiBB, 
the voltage drop ratio of three barriers, {i = 1,2,3), 
is required. Figure 5(a) shows dlso/dYsD at 12T and 
we mark L1-L4 for the current threshold lines from the 
N = 1 Coulomb diamond. At the threshold LI, the Fermi 
energy of the left electrode is aligned with the energy level 
of the left dot. The current threshold several orders of 
magnitude smaller, marked by L2, indicates the onset of 
a cotunneling process where the right-dot energy level is 
aligned with the Fermi energy of the left electrode. Simi- 
larly, the Fermi level of the right electrode is aligned with 
the right (left) dot at the threshold line marked as line 
L3 (L4). The slopes of these threshold lines L1-L4, that 
is, -32, -54, 46 and 62, respectively, give (i = 1, 2, 3) of 
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FIG. 5: (a) dlso/dYsD plotted in logarithmic color scale 
in \ow-VsD region measured at 12 T. Four current threshold 
lines of the N = 1 Coulomb diamond are marked as L1-L4. 
The slopes of these lines gives the voltage drop ratio of the 
triple barrier, ai (i = 1,2,3), for the \ow-Vsd region as 0.27, 
0.19, and 0.54. (b) Resonant current peak under 39.3 GHz 
microwave application to Vsd- Interdot microwave- assisted 
tunneling is clearly apparent. Measurements were performed 
at 1.5 K using a different sample with the same triple bar- 
rier structure wafer. Distance between the main and satellite 
peaks was used to estimate the interdot voltage drop ratio 
a2 in high-l/^i:) regions as 0.193. The microwave powers are 
labeled for each trace. Each trace is shifted by a constant 
value. 

0.27, 0.19 and 0.54, respectively, at Vsd -^0. Figure 5(b) 
shows the effect of microwave irradiation on the resonant 
current peak. With increasing microwave power, satel- 
lite current peaks appear on the left (right) side of the 
main peak, owing to interdot tunneling upon the absorp- 
tion (emission) of microwaves. Such microwave-assisted 
tunneling in a double dot has been observed in lateral 
and vertical dots [l^, [l^. At the satellite peak posi- 
tion, interdot detuning is equal to the microwave photon 
energy. Thus the distance between the main and the 
satellite peak (marked as two dotted lines in Fig. 5 (b)) 
gives a2 of 0.193, which is close to the value ioiVsD 0. 
Similar values are obtained for different current peaks at 
different Vsd''^ from up to 100 mV, indicating that 0^2 
is independent of Vsd- On the other hand, voltage drop 
ratios for outer barriers {i = 1,3) are found to depend 
on VsD' Indeed, in Fig. 1, the current threshold line from 
= at positive Vsd is not straight but bends upward. 
This indicates decreasing ai with increasing Vsd- As the 
slope of the current threshold lines is -12, as determined 
from Figs. 2(a)-(d), ai at Vsd ~35 mV is smaller than 
ai at Vsd by 12/32. Thus {i = 1,2,3) is 0.10, 
0.19 and 0.71, respectively, at Vsd ^^35 mV [11]. Us- 
ing these a^'s as well as the slope dit Vsd ^^35 mV, we 
obtain of 0.89 for our Ino.04Gao.96As dot, and the 
g-factor difference between the two dots is 0.56. Thus, 
l^i^l for our GaAs dot is 0.33. The value for the GaAs 
dot is consistent with that obtained in the previous work 



on a vertical single dot with the same 10-nm-thick GaAs 
well p^]. The value of 0.89 for InGaAs seems to be large 
compared with the value obtained in the study of the g- 
factors of bulk In^^Gai-a^As [20], probably owing to the 
increasing density of In at the center of the well rather 
than at the InGaAs/ AlGaAs interfaces, as is seen in the 
self- assembled InGaAs dot systems [21]. A stress in the 
thin InGaAs layer can be another possible reason for large 
g-factor. 

The width of the current peak line in both regions A 
and D is ~0.5 mV in Vsd^ and corresponds to ~0.1 meV 
using 0^2 • This width gives an upper bound for possi- 
ble contributions from phonon- absorbing/emitting tun- 
neling, and is smaller than the Zeeman energy difference 
in high magnetic field (^10 T). Thus phonon absorp- 
tion/emission for a misaligned spin sublevel in region B 
(C) in Fig. 3 is negligible. This point is also confirmed 
by a more detailed calculation, as in Fig. 3, where the 
phonon effects are included[16]. As depicted in the right 
inset of Fig. 1, there are some unoccupied levels located 
a few meV's below the aligned level in the right dot. 
Thus, in addition to the resonant tunneling through the 
two aligned levels, there is always a tunneling process 
with meV-phonon emission down to these low-lying lev- 
els. Such meV-phonon-emitting tunneling gives a back- 
ground current of the order of ~10 pA, which defines the 
current threshold from the N = region, as seen in Fig. 
1 and Figs. 2(a)-(d). Resonant tunneling current peak 
lines are also seen for negative Vsd- However, most of the 
peak widths are large and no clear kink structure is ob- 
served. For negative Vsd^ the direction of the A/2 shift 
would be opposite, and Si of the kink should be smaller 
as it is governed by the g-factor in the right dot, which 
is smaller for GaAs. This may be one of the reasons why 
we were not able to detect clear kink structures. In Fig. 
2(d), there seems to be an additional line parallel to the 
current threshold line from the N = region. We found 
that these lines, which appeared irregularly with some 
magnetic fields and Vs-d's, may be due to a fluctuation 
of the density of states in the electrode [22^. 

In conclusion, we measured the single electron trans- 
port through two zero-dimensional states with different 
g-factors. We found that the resonant tunneling is sup- 
pressed even when one of the Zeeman sublevels is aligned 
within the transport window. Finite level broadening of 
zero-dimensional states partially relieves this bottleneck 
effect, and gives a current peak when level detuning is 
set to half the Zeeman energy difference. 
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